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Abstract 

The front form Hamiltonian for quantum chromodynamics, reduced to an effective Hamil- 
tonian acting only in the qq space, is solved approximately. After coordinate transformation to 
usual momentum space and Fourier transformation to configuration space a second order differ- 
ential equation is derived. This retarded Schrodinger equation is solved by variational methods 
and semi-analytical expressions for the masses of all 30 pseudoscalar and vector mesons are 
derived. In view of the direct relation to quantum chromdynamics without free parameter, the 
agreement with experiment is remarkable, but the approximation scheme is not adequate for the 
mesons with one up or down quark. The crucial point is the use of a running coupling constant 
a s {Q 2 ), in a manner similar but not equal to the one of Richardson in the equal usual-time 
quantization. Its value is fixed at the Z mass and the 5 flavor quark masses are determined by 
a fit to the vector meson quarkonia. 
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1 Introduction and motivation 



One of the most outstanding tasks in strong interaction physics is to calculate the spectrum and the 
wavefunctions of physical hadrons from quantum chromo dynamics (QCD). Discretized light-cone 
quantization (DLCQ) pj has precisely this goal. Its three major aspects are: (1) a rejuvenation 
of the Hamiltonian approach, (2) a denumerable Hilbert space of plane waves, and (3) Dirac's 
front form of Hamiltonian Dynamics. In the front form Q, or in light-cone quantization ||, one 
quantizes at equal 'light-cone time' x + = t + z, as opposed to the conventional instant form where 
one quantizes at equal usual time t. As reviewed in the front form has unique features, among 
them: The vacuum is simple, or at least simplier than in the instant form, and the relativistic 
wavefunctions transform trivially under certain boosts [& 0]. Both are in stark contrast to the 
conventional instant form. Over the years, the light-cone approach || has made much progress. 
Calculations [||] agree with other methods particularly lattice gauge theory. Zero modes of the 
fields can be important carriers of quantum structures Q, particularly of those of the vacuum |7[ S|. 
Dimensionally reduced models @, || provide much insight into the structure of possible solutions 
to QCD. But chiral aspects are not yet understood, and non-perturbative renormalization remains 
a challenge || |l(]] as for any Hamiltonian approach. 

But despite the many successes of light-cone Hamiltonian methods one misses the contact to 
phenomenology beyond the perturbative regime. We believe that more QCD-inspired approaches 
are needed, work like for example [11, 12 1 or [[b^, 14], where the formalism is related to the experi- 
ment. The present work is of this type. 



Right from the outset when applying DLCQ to gauge theory in 3+1 dimensions [15, 16, 17] 



it was clear that one should need an effective Hamiltonian. In [|R| an integral equation in the 
light-cone momenta was solved numerically, which was derived by procedures similar to those of 
Tamm and Dancoff and a non-integrable singularity was removed by an ad hoc assumption. 

But recently |^] the method of effective interactions was generalized to avoid the usual truncation 
in the particle number |21]. As it turns out, one can assemble all many-body aspects into a vertex 
function which bears great similarity with the running coupling constant. 

One wonders: How can such a simple structure account for the spectra and wavefunctions of 
all scalar and vector mesons? Is this not too much of a claim? On the other hand the effective 
Hamiltonian has been derived [21] from the QCD-Lagrangian without condition on the coupling 
constant or on the mass of the constituents. One way of checking this is to compare to experiment, 
and this shall be done in this work very roughly and preliminarily. Lacking the running coupling 
constant going with the theory pl|, one can replace it by one of its current phenomenological 



versions ]2^, g4j. The present work applies the one of Richardson [22]. It interpolates smoothly 



between asymptotic freedom [25, Efl] and infrared slavery. After that, one has no freedom in 



the theory and no adjustable parameters. Since the quark masses cannot be determined from 
independent measurements, they must be determined self-consistently from a fit to some of the 
meson masses. This in itself is not trivial, except when having analytical expressions. 

In particular, a coordinate transformation from front to instant-form coordinates is performed 
in Section 3. Apart from a more transparent interpretation, this way of writing down the integral 
equation has certain advantages in performing the calculations. No assumptions will be made in this 
section: All manipulations are straightforward and fully equivalent to the front-form formulation. 
In Section 4, the bound-state equation is approximated semi-relativistically which allows for Fourier 
transforming the momentum-space integral equation into a configuration-space Schrodinger type 
equation. The so obtained Hamiltonian is reduced in Section 5 to a minimal number of terms 
(Coulomb plus linear potential plus one spin-dependent term distinguishing between singlet and 
triplet) and diagonalized approximately by a variational method. 

The masses of all pseudoscalar and vector mesons in Section 6 are thus semi-analytic and 
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Figure 2: The quark-antiquark potential V(r) 



Figure 1: The effective interaction in the qq sec- 
tor. By exchanging 'effective gluons', a single 
quark state with four-momentum ki and spin 
projection Ai is scattered into the quark state 
(k[, X[). Correspondingly, the antiquark is scat- 
tered from (k2,\2) to (k^X^)- 



as given in Eq.(ll) is plotted versus the relative 
distance r. The Coulomb and the confining po- 
tential are indicated. Some of the heavier meson 
masses are inserted to provide a scale. 



approximate solutions to a second order differential equation in configuration space. In comparison 
with the empirical masses [p8f , they are not much worse than those of potential models 30, 31, 
32, 33 1, or predictions from heavy quark symmetry ]3^| , or even predictions based on lattice gauge 
calculations [35, 36, 37]. In view of the direct link to QCD [^] and the simplicity this should be 
regarded as considerable progress in the front-form approach. 

But there is a potential danger in such an endeavor. The present work is motivated by the 
question whether the simple structures to be displayed can describe experiments at all. Obviously 
they can, but it should be emphasized that numerically accurate solutions need another effort. This 
is currently being attempted [38] and has a different objective than to develop models designed to 
reproduce the data. 



2 The effective Hamiltonian for QCD 

In Discretized Light-Cone Quantization (DLCQ) one seeks to solve the eigenvalue problem 

#L C |*&) = M^ b ) (1) 

for a field theory. The 'light-cone Hamiltonian' -£Tlc = P^Pa @ is the Lorentz invariant contraction 
of the energy-momentum four-vector and has the dimension <mass> 2 . The eigenvalues M 2 are 
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interpreted as the square of the invariant mass of state b. Working in momentum representation, 
the three space-like components of P 11 are diagonal operators, with eigenvalues P + = J2j and 
P± = J2j k-Lj ■ The sum runs over all particles in a Fock state. Each particle has a four-momentum 
denoted by kj = (kj ,k± j: kj) and sits on its mass shell (k^k^j = m 2 . The time-like component, 
the Hamiltonian proper P~ , is a complicated and off-diagonal operator acting in Fock space. Its 
matrix elements are tabulated in ||. Based on the boost-properties of light-cone operators [|| 
one can transform to a frame where P± = 0, thus P^P^ = P + P~ ■ Since P + is diagonal, the 
diagonalization of P~ and of i?LC amounts to the same. The Hilbert space for diagonalizing P~ 
is spanned by all Fock states which have given eigenvalues of P + and P± = (0, 0) and can be 
arranged into sectors according to the particle number like qq, qqg, or qqqq. For any fixed value of 
the harmonic resolution K = 2LP + /it, the Hamiltonian matrix in Eq.(jl]) is finite and in principle 



could be diagonalized numerically [y. Details can be found in the literature 20, 21]. 

DLCQ is quite useful to cleanly phrase the problem, but to do calculations particularly in 3+1 
dimensions one has to develop effective Hamiltonians. Fock space truncation in conjunction with 
perturbation theory in the manner of Tamm and Dancoff [18, 19] is unsatisfactory, because one has 
to resort to ad hoc prescriptions to make things work fl~6[| . These drawbacks can be avoided by the 
method of iterated resolvents p0| , ^l| ] . It turns out possible to convert the many-body matrix Eq. (|l]) 
into a well defined two-body equation with an effective interaction acting only in the qq-sp&ce, i.e. 
Hcfi\ipb} = Mfr\t/jb). In the continuum limit one has to solve the integral equation 

M 6 2 (x,4;Ai,A 2 |V' 6 )= £ /^ , d 2 ^(x ! ^;Ai,A2|^effk , ,fcl;Ai,A^)(x / ,fc|;Ai,A^|^). (2) 

A 1 ,A 2 

The bras and kets refer to qq Fock states which can be made invariant under SU (iV) , 

1 n c 

|x,4;Ai,A 2 ) = Ai)4(Afe,A 2 )|0). (3) 

Goal of the calculation are the momentum-space wavefunctions {x, k±; Ai, A 2 |V ; b)- They are the 
probability amplitudes for finding the quark with helicity projection Ai, longitudinal momentum 
fraction x = kf /P + and transversal momentum k±, and correspondingly the antiquark with A 2 , 
1 — x and — k±. The effective interaction as diagrammatically displayed in Figure [l] is a sum of three 
terms: The first two diagrams are kind of a one-gluon exchange and describe the flavor conserving 
part of the effective interaction, while the last graph due to the two-gluon annihilation can change 
the flavor. In the present work we deal only with the first of them. The kernel of the integral 
equation (|2|) has a diagonal 'kinetic' and an off-diagonal 'interaction' energy, i.e. 



Mi (x,k x ;Xi,X 2 \^b) 



mi + kj_ m 2 + kj_ 
x 1 — x 



(x,k±;X 1 ,X 2 \ipb) 



VK l ,j-f,J Q z vx(l - x)x'(l - x') 

A 1> A 2 

The most important factors are the four-momentum transfer 

Q 2 = ~(h - k[) 2 = -(k 2 - k' 2 f (5) 
and the vertex function r(Q, A) which likes to combine with the coupling constant g to become 
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the like-to-be 'running coupling constant'. For QED, this factor reduces to the fine structure 
constant (3 = a ~ 1/137. The spinor factor S(Q) represents the familiar current-current coupling 



(A 1 ,A 2 |5(Q)|Ai,A' 2 ) = [u(fc 1 ,Ai)y 4 w(fc / 1 ,A , 1 )] [u(k 2 , A 2 ) 7m u(A; 2 , A 2 )] . 



(7) 



The cut-off function Q(x',k\ ) restricts integration in line with Lepage-Brodsky regularization 



@(x, k± 



m\ + k[ 



+ 



rn 



2 + k 2 

2'± 



X 



1 



X 



< (mi + m 2 ) 2 + A 2 



(8) 



The mass scale A can be chosen freely. 

Despite having been derived in the light-cone gauge A + = 0, the effective interaction is man- 
ifestly gauge invariant, depending only on the quark currents. The instantaneous interaction has 



cancelled exactly against other gauge-variant terms, see |15, 16, 21]. Since one works in the front 



form, it is also frame and boost invariant. Explicit calculations for QED [16, 38] are numerically 
very stable, and reproduce quantitatively the Bohr spectra and the fine and hyperfine structure. 
The vertex function hidden in the like-to-be running coupling constant of Eq.(||) has the same 



perturbative series expansion as the running coupling constant [21] which is indicated in an artists 



way in Figure |l|. What is missing, thus far, is a renormalization group analysis of the formal 
expressions. In the absence of that, we are interested in consequences of Eq.(|5|). How can it be that 
such a simple expression accounts for hadronic phenomena? What are the invariant masses of the 
pseudoscalar and vector mesons, using such an interaction? How far does one get with analytical 
procedures, and in particular where does the approach go wrong? 



Lacking an exact expression for a s (Q 2 ), one can resort to reasonable parametrizations 22, 23, 
24]. In the sequel, we shall content ourselves with the form of Richardson 22], 
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27 \n{a 2 + Q 2 /n 2 ) 



(9) 



At least, this form interpolates smoothly between asymptotic freedom 25, |2j|] and infrared slavery. 
In the original work, the the parameter a was set to have the value a = 1 and n was kept as free 
parameter to be determined by the spectra. Here, we take the value of a s (Mz) = 0.1134 ± 0.0035 
as measured at the Z-mass to fix 

K = 193MeV. (10) 



Its Fourier transform 22] generates two terms, see also below, 



8tt / 1 
27 V r 



+ k 2 r 



(11) 



a (strong) Coulomb and a linearly rising potential, as plotted in Figure versus r = \x\. The 
linearity of the confining potential is a consequence of a = 1, as used in p2| and throughout the 
present work. If one varies a one gets the curves displayed in Figure |3|. It is taken from [27]. Here, 
we do not want to keep a as a free parameter. For one reason, we refuse to speculate at this point 
whether the potential is strictly confining or not. For the other reason, the results to be displayed 
below are not very sensitive to large distances since the wavefunctions decay rapidly. Last but not 
least, one has to await the renomalization group analysis of the a s (Q) which truly comes with the 
theory @. 

The flavor quark masses are then the only free parameters of the approach. Of course, they 
are subject to be determined consistently by experiment. Natural candidates are the masses of the 
pseudoscalar (0~) and vector mesons (1~). Since the flavor quark masses potentially range from 
a few MeV up to some 100 GeV, see for instance Figure 0, one faces two problems: (1) By good 
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Figure 3: The confining potential as function of 
the parameter a. Values are from top to bottom: 
a = I, a = 1.0005, a = 1.01, and a = y/e ~ 1.65. 
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Figure 4: Bound states of qq-pairs versus quark 
masses. All masses are given in units of the 
QCD-scale k. The upper curve refers to the 
triplet (S e = 2), the lower to the singlet 
(S e = 0). The masses of some vector mesons 
(p°, co, 4>,J/ip) are marked by a (A ), those of the 
some pseudoscalars (tt°, r], rf ,rj c ) by a (*). 



reasons, the numerical solutions of the integral equation have been restricted thus far to systems 



with equal masses of the constituents like positronium [16, 38 1. The wave function is then peaked 
at x = m\/{m\ + m-z) = 1/2. For very asymmetric systems this will be a problem. To avoid 
that, we shall identically rewrite the integral equation in the next section in terms of instant-form 
variables, which are somewhat easier to deal with. (2) Shall one really perform a calculation of 
similar complexity as in preceeding work [|l^, ^] for any given set of m\ and m,2 when intending 
to fit them to the meson masses, or shall one aim for a quasi-analytic but approximate solution? 
In view of the preliminary character of the present study, we have opted for the second. At least, 
this will pave the way for a future and improved solution. 
In the sequel, we shall replace H e s by the operator 

H= %^T^)( H '-^ + m ^) ■ (12) 

It differs from H e Q by an additive constant and an overall scale, which both are Lorentz scalars. 
Both H and its eigenvalue E have the dimension of a (mass) and have much in common with 
the non-relativistic Hamiltonian and the binding energy E, as we shall see. As compared to an 
instant-form Hamiltonian, however, the main advantages of the front-form H e g as given in Eq.(j5|) 
prevail, namely the additivity of the interaction and the Lorentz invariance of the eigenvalues. - 
The rest of this work is a simple and straightforward evaluation. 
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3 Transforming variables from the front to the instant form 

The single-particle four- momenta can be parametrized either in the front form, k± = (kf ,k±, k7) 
and k 2 = {k 2 , —k±, k 2 ), or in the instant form, 

K = (k%, k±, k z ) = {Et,k), and k% = (k° 2 , -k ± , -k z ) = (E 2 , -k) . (13) 
Since k^ki^ = mf, the time-like components are functions of the space-like components 

k - = m i+ k * } or Ek = yjm 2 + + k 2 = \j m 2 + k 2 . (14) 

The transformation function between x and k z is obtained straightforwardly from P + , i.e. 

x = x{k z ) = |jr^jr • ( 15 ) 

The front form integral equation is boost and frame invariant, and therefore can be solved also in 
the center of mass frame, where the total momentum P vanishes. Changing integration variables, 
Eqs.@ in conjunction with Eq.(|T^) can thus be rewritten identically as 

E (A#) = r(jfe) (£#) + [ d 3 k' (k\U\k') (k'\i/>). (16) 



For simplicity the explicit summation over the helicities is suppressed. Contrary to Eq.(||) all three 
integration variables have now the same support. The kinetic energy 

T{U) = 2(m 1+ m 2 ) l,^^ + ~T^T ~ ) ~ 2(m 1+ m 2 ) 

becomes the familiar expression with the reduced mass m r , for sufficiently small momenta, 

T(k) = ^- , with — = — + — . (18) 
2m r m r mi ^2 

Nevertheless, there are explicit residues from the front form. The Fock state \x, k±) has the same 
P + as \x', kj_). Expressed in instant form variables P + = P° + P s . Since P 3 = k z i + k Zj2 = one 
is left with P° = P ' , or explicitly 

Ei + E 2 = E[ + E' 2 , or k 2 = k' 2 (19) 



for every matrix element. Obviously, the interaction kernel in Eq.fllq) cannot change the size of k, it 
only changes its direction. This is a source of great simplification. For example, the four-momentum 
transfer is always identical with the three-momentum transfer 

Q2 = _ {h _ k r r{h _ k r )fM = {ki _ £,)2 _ {Ei _ E r )2 = ^ (2Q) 
and the three- momentum transfer and its mean, 

q = k-k' and p = -(k + k'), (21) 
respectively, are always orthogonal: 

p. q=\{k + k')(k - k') = \{k 2 - k' 2 ) = 0. (22) 
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The Jacobian of the transformation Eq.(|T5|) is evaluated by means of the identities 

dx _ (E 1 + k z )(E 2 - k z ) ( El + kz) (E 2 -k z ) 

dk z ~ E 1 E 2 (E 1 + E 2 ) and x ^- x >- (E 1+ E 2 y • (23) 

The auxiliary functions 



A(k,k') = J )Z 1 ' ?( , and B{k,k') = (^ + ^) (24) 



are useful for factorizing 



' (E 1 + k' z )(E 2 - 


K) 


(E 1 + k z )(E 2 - 


k z y 


m r dx' 




\Jx(l — x) x'(l 


- x') 



E\ E\ 



A(k,k')B(k,k') dk z , (25) 



and to single out A which is not rotationally invariant. Both A and B are dimensionless and of 
order unity for sufficiently small momenta. The cut-off function G', as introduced in Eq.(||) to 
define a maximum transversal momentum restricts of course also three-momentum: 

y J ~ \2 J A 2 + (mi + m 2 ) 2 v ; 

The quark currents in Eqs.(||) or ([/]) can be evaluated with the Gordon decomposition []39|]. Since 
we work in the Lepage-Brodsky convention for the spinors § one has 

u(fci, Ai)7u(fc / 1 , A' x ) = (2p - i q x o\) AiA , , (28) 
with .R = iqx p. (29) 



These expressions are simpler than usual, due to Eq,(|l9|). For the antiquark one must change the 
sign of both k and k' , and replace the quark-spin matrix a by r. The current term becomes then 
explicitly 

J(h,h') = j^^- Ai) 7 " »M>A'i)] WhMwi&X)] (30) 

- c ?™ g-^. + ,„ Ui + g- f2/ l + *- f ^ 



(£i+7rii) 2 (^2+m 2 ) 2 / \ (.Bi+mi) 2 (E 2 + m 2 ) 2 J 
_ ( - 2f- iqx f ), (31) 



i?i + mi / V E 2 + m 2 

with C = (^+^)(^ + ^). (32) 
Amim 2 

After A and B a third auxiliary function C is introduced, which is also dimensionless and of order 
unity. As expected for a Lorentz scalar, J is rotationally invariant. 

Thus far, all quantities considered are of order unity for sufficiently small momenta. The most 
important part of the interaction kernel 

(k\U\k') = e(fc') A(k, k') B(k, k') J{k, k') V(k, k') (33) 

is therefore the interaction proper 
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It depends only on the three-momentum transfer. 

The front form is frame and boost invariant, as mentioned. It is rotationally co- but not 
rotationally invariant, particularly when the spatial rotations are performed perpendicular to the 
z-axis. This aspect is reflected in the appearance of the factor A as defined in Eq.(24). The 



violation of rotational invariance occurs however in such a form that it can be absorbed into the 
wavefunction. If one inserts 



into Eq.(16), the factor A cancels in the new integral equation 

E <p(k) = T(k) 4>(k) + f d 3 k' Q(k') B(k, k') J{k, k') V(k, k') </>(&) . (36) 



The kernel is now rotationally invariant. Since no approximations have been made, the solutions 
of this equation, mutatis mutandis, are identical with those obtained from the original front-form 
integral equation, Eq.(||), but Eq.(|36[) is much easier to deal with. 

4 The retarded Schrodinger equation 

The front form of Hamiltonian dynamics [|| has wonderful properties but it does not appeal strongly 
to our intuition, not even when it is transcribed to instant form variables. Thinking in terms of 
momentum-space integral equations is not always easy. The equations become more transparent 
when Fourier transforming them to configuration space and the corresponding Schrodinger form of 
quantum mechanics. 

We begin with rewriting Eq,(|3"6|) conveniently as 

E(f)(k) = Jd 3 k'H(q,p) <j)(k') . (37) 

The kernel H is expressed in terms of the momentum transfer and its mean rather than by k and 
k' . It is the Fourier transform of the Schrodinger Hamiltonian. To see that one multiplies the whole 
equation with exp(ik • x) and integrates over d?k. Defining the Fourier transforms by 

= Jd 3 k 4>(k') , and H(x,p)= Jd 3 qe ic i g H(q,p) , (38) 

one gets an eigenvalue equation of the Schrodinger type with a possibly non-local Hamiltonian 

Eip(x) = H(x,p) ij)(x) , with p= —iV x . (39) 

The momentum transfer q is Fourier conjugate to the position x of the quark in the center-of-mass 
frame, and p is the associated momentum operator. This holds in general, but unfortunately one 
is unable to perform the Fourier transform explicitly with all the square roots behind the energies 
Ei. The way out is, of course, to expand and to develop a systematic approximation scheme. We 
base it on the Lepage-Brodsky cut-off and choose A such that 

^ < 1 , (40) 
for the lighter quark m\. All square-roots are expanded to first non-trivial order 

Ei oi mi + = mi + + , (41) 
2m,- 2m,- arm 



9 



which is a semi-relativistic approximation. In the worst case it allows for relativistic velocities of 
the lighter particle up to \k\ ~ mi. The expansion of the various factors in the kernel of Eq.(|36|) 
yields up to second order 

C = l + + with +-L (43) 



f> 2 
2m 2 


Q 2 

8m 2 ' 


with 


1 

m 2 


p 2 
Am 2 a 


q 2 
16ml ' 


with 


1 

ml 


P 2 


Q 2 


(a x 


q) ■ (r 


1m\m2 


8m 2 




Am\m2 



S-R + t-R S-R 

Am 2 Am 2 , m\m2 



respectively. One should emphasize that the form of V(q) need not be known at this point, since 
it does not depend on p. The total spin and the kinetic energy, 

S + m d (45) 

respectively, complete the definitions. Finally, one can conjecture that the wavefunction decays 
sufficiently fast, such that it acts itself like a cut-off. We therefore set Q(k') = 1. 

The Hamiltonian operator in Schrodinger representation becomes then straightforwardly 

H - JL ( ! + J4U t + V(r) + + <" f • g X W) 

2m r V mi + m2/ - 8mj Am\m2 

18V fa-L t ■ L S • L \ 



r 



with the usual angular momentum operator L = x x p. Since the average potential is spherically 
symmetric, one uses (cr x V) • (f x VV) = §(<?• t)V 2 V" and a ■ f = 2S 2 — 3 to get 

H = 1 ( 1 | ^ ) f | r (r) 3m r (m 1 -m 2 )^ 



2m r V mi+m2/ — 8m 2 \mi+m2 (m\ + m2) 2 ) 

ldVfa-L t • L S-L\ V 2 V ~ 2 
+ —5- — -» + -;— r + o 5 • 47 

r or \ 4m f Am% m\m2 I 6m\m2 

Its structure is a direct consequence of gauge theory particularly QCD and holds for an arbitrary 
running coupling constant a s (Q 2 ). We emphasize particularly that this structure was obtained 
from a fully covariant theory The statement could even be stronger without our inability to 
evaluate the Fourier transforms without expansions. Richardson's parametrization of a s (Q 2 ) yields 
the potential V(x) as given in Eq . (|TT|) , and thus 

V 2 V(r) = p(—+A^S(x)) and = ( - + \) , (48) 

V r I r or \ r r° I 

with = |f ~ 0.93. If one works with QED, one sets k = and chooses the value ~ 1/137. 

We now have reached our goal: The retarded Schrodinger equation and its Hamiltonian have a 
wonderfully simple structure which can be interpreted with ease. The average potential V(x) plays 
a different role in the different terms of the equation. In the first term of Eq.(47), in the kinetic 



energy, it generates an effective mass of the quark which depends on the relative position and which 
reflects the non-locality of the interaction. In the second term, V(x) appears in its natural role 
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as a potential energy. In the third term one observes V{x) as the analogue of the Darwin term. 
In the remainder V(x) provides the coupling strength for the analogue of the fine and hypefine 
interactions of atomic physics particularly the spin-orbit interaction. Contrary to common belief 
they exist not due to weak coupling but they appear also for strongly coupled QCD. 

Finally, one must come back to the expansion scheme of Eq.(41). Its validity cannot be jugded 



a priori, since the expansion is made under the integral. The omitted terms are of second order in 
p 2 /m 2 for the lighter quark. Whether this is justified or not can be decided only a posteriori, by 
the expectation value of the omitted next higher term 



i f(p 2 )\ 2 



2 

mf 



(49) 



Only if 5 is (very) small as compared to unity, the expansion in Eq.(41) is justified. If it is 
comparable or larger than 1, the solution must be rejected, and another regime of approximation 
must be found. Below, we shall meet cases like that. 



5 Meson masses by parametric variation 

It will take some time and effort to work out all the many consequences of the integral equations, 



Eq.(|j) or Eq.(pq), or of the retarded Schrodinger equation (47). In the sequel, we shall restrict 
ourselves to calculate only the ground-state masses of the pseudoscalar and vector mesons. If one 
leaves aside the recently discovered top quark and restricts to 5 flavors (u, d; s, c; b), one has thus 
30 different physical mesons, since charge-conjugate hadrons have the same mass. 

One cannot calculate these masses, however, without knowing the quark mass parameters mi 
and 77i2. These cannot be measured in a model-independent way. In the sequel, we shall adopt the 
point of view that they have to be determined consistently within each model, for the better or the 
worse. One has thus 5 mass parameters to account for 30 physical masses. Which ones should be 
selected to fit? There are 142 506 different possibilities to select 5 members from a set of 30, and 
we have to make a choice: We choose the five pure gg-pairs. Even that is not unique: Shall one 
take the pseudoscalar or the vector mesons? We shall do both. 

Of course, one runs into the problem of the chiral composition of the physical hadrons. In order 
to avoid that in the very crude estimate below we shall substitute all physical mesons by pure 
gg-pairs - by fiat. We shall thus replace the 'pions', for example, by 'quasi-pions' with the same 
physical mass. The uu-, ud-, du-, or dJ-eigenstates shall be identified with the quasi-7r°, quasi-7r + , 
quasi- 7r", or the quasi-77, and so on. This simplification will be revoked in future work and is, by 
no means, a compelling part of the model. 

Our problems lie in another ball park. One should not deal head-on with the full complexity 
of the integral equations, or of the retarded Schrodinger equation. Which part of the Hamiltonian 
should one select in a first assault? Some help is gained by the rather unique property of the light- 
cone Hamiltonian: kinetic and interaction energy are additive. One can select always an 'interesting 
part' Ho, and check e posteriori, by calculating the expectation value of AH = H — Hq, whether 
the selection makes sense, or not. Since the scalar and pseudoscalar mesons have only little orbital 
excitations, i.e. are primarily s-waves, one can disregard the spin-orbit part and choose first 

1 V(r) \ p 2 V 2 V ^ 2 V 2 ^ ( 3m r _ (mi - m 2 f ' 



go- 1+ — H— + y(r) + -^—S 2 -^[ - 1,1 . (50) 

\ m\ +7712/ 277i r 3777,17772 OT77^ \ 777i + 771-2 \ m l + ^2} J 

Even that looks to complicated for a start-up. We therefore select those terms which have turned 
out to be important in the past, namely the central potential and the triplet-interaction mediated 
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Table 1: The flavor quark masses in MeV, as obtained from a fit to Eq.(57). The first row refers 
to a fit for the singlets, the second to the one for the triplet. 



Flavor mass 


u 


d 


s 


c 


b 


From fit to 0~ 


2.3 


155.6 


430.6 


1642.3 


5330.8 


From fit to 1~ 


222.8 


236.2 


427.2 


1701.3 


5328.2 



by the total spin. Omitting the effective mass and the Darwin term, our choice is therefore 

P 2 rrt \ 2K ^ P 2 I 3 2k2 S 2 o 2 

H Q = ^— + V{r) + = - tL + + p K \ . 51 

Working in a spinorial representation which diagonalizes S z and S 2 , we replace the latter by the 
eigenvalue S e = S(S + 1) and take or 2 for the singlet or the triplet, respectively. 

How does the wavefunction for the lowest state look like? For a pure Coulomb potential the 
solution has the form 

V(x) = -^Aie" Ar . (52) 

\/7T 



Omitting the Coulomb part, a linear potential can be solved in terms of Airy functions and its 
integral transforms [40]. If one has both, one will have some mixture of the two. But for the 
present start-up check even that requires too much effort. 

Rather shall we pursue a variational approach and choose Eq.(|5l]) as a one-parameter family 
(A). One could take also harmonic oscillator states fi"3| , |i"4|| , but with Eq.(|5l"l) the expectation values 
are particularly simple: 

(V| f |V) = A 2 , (VI - IV) = A , and (V| r |V> = |r • (53) 

r 2A 

These are all one needs for calculating the expectation value of the energy 

E = (VlilolV) = f- ~ /3A + ^ + ^f\- (54) 

Since we deal only with ground states, we are not in conflict with the statement that the wave- 
function cannot be purely Coulombic. For the pure Coulomb case, the IS- and IP-states would 
be degenerate and the respective ratio |V2s(0)| 2 /|Vis(0)| 2 = 0.125 would disagree with the experi- 
mental values ~ 0.63 and ~ 0.50 for charmonium and bottomium [pS[|. 

We aim at calculating the total invariant mass of the hadrons and return to the light-cone 
Hamiltonian Hlc = M 2 , i.e. to M 2 = (mi + 1712) 2 + 2(mi + m2)E. For equal masses mi = 771,2 = m 
one preferably expresses the variational equation ( ]54| ) in units of the fixed QCD scale k, introducing 
the dimensionless variables 

A ^ 777 , „ T fM\ 2 



£ = - , and W=[ — \ . (55) 



The variational equation (p3) reduces then simply to 



in>:() ^-(^- 2 -jf)s + Z 2 + ^- s - (56) 
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Table 2: The validity check. The first row displays the values of (A/m) 2 as obtained from the mass 
fit to the singlets, the second row those from the mass fit to the triplets. If the mean momentum 
is comparable to the mass, or larger, the solution has to be rejected. The extremely large value for 
the u-quark in pseudoscalar fit gives a good example for such a case. 



A 2 
m? 


u 


d 


s 


c 


b 


For fit to 0" 


293 


1.45 


0.53 


0.25 


0.22 


For fit to 1" 


0.79 


0.76 


0.48 


0.25 


0.22 



We must vary A, thus s, such that the energy is stationary, 



dW(s;0 



ds 



*(0 



0, thus = W(s*(0) = W*(0 , (57) 



at fixed values of the parameters (£, /3, S e ). This leads to a cubic equation in s which can be solved 
analytically in terms of Cardano's formula. In special cases they can well be approximated by a 
quadratic equation, namely when s* 3> 1 or when s*<l. We got accustomed to refer to these two 
regimes as the Bohr and the string regime, respectively. In the Bohr regime the Coulomb potential 
dominates the solution and the linear string potential provides a correction. In the string regime 
the linear string potential dominates, with the Coulomb potential giving a correction. Solutions in 
the string regime, however, imply that the ratio (p 2 )/m 2 = A 2 /m 2 becomes so large that one is in 
conflict with the validity condition Eq.([49|). 

Rather than to display explicitly the straightforward but cumbersome formalism, we present the 
analytical results in the graphical form of Figure |j. The total mass M = 2k\/W* is almost linear 
in the quark mass, with small but significant deviations. In line with expectation, the hyperfine 
splitting decreases with increasing quark mass. Less expected was that the splitting increases so 
strongly with decreasing quark mass. For very small quark masses, the triplet mass starts off 
at a finite and almost constant value. The singlet mass takes off from zero like a square root but 
unfortunately not linearly as required by the soft pion theorems. Determining the u-mass by fitting 
to the quasi-pion gives a value close to the 'current mass', see Table |l|. The resulting s<l, see 
Tables |l| and ||, implies the ultra-relativistic string regime and that the validity condition is badly 
violated. The retarded Schrodinger equation with its semi-relativistic approximation scheme is thus 
not appropriate for describing quasi-pions. For the n and the rj' , the scaling variable s is of order 
unity, while for the quasi-ry c one definitely is in the Bohr regime. Here the masses are similar or 
close to what is refered to as the 'constituent-quark' mass. 

Since singlet and triplet are so close for s > 1, one fits the quark masses preferentially with 
the vector mesons. In the lack of empirical data we have set M„, = Mx, which should be of minor 
importance in the present model, see Figure |I[ The flavor masses are now in close agreement 
with the constituent quark masses, see Table [l], and the smallness condition is satisfied better, see 
Table |. 

Having determined the quark masses one has exhausted all freedom in the model. We now ask: 
How well agree the remaining 25 meson masses with experiment? - The formal procedure runs 
quite analogously, except that it is now easier. With the masses fixed, one varies A separately for 
each flavor composition subject to Eq.(^4f). The results are compiled in Table and compared 
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Table 3: The masses of qq-hadrons are compared with experimental values. The flavor quark 
masses used are inserted in column 2 and come from a fit to the vector mesons. The first line 
within each box refers to the hadronic symbol of the meson; the second line gives the calculated 
(measured) vector mass in MeV; the third line accounts for the calculated (measured) pseudoscalar 
mass in MeV, and finally the fourth line specifies the hadronic symbol of the pseudoscalar meson. 





m q 


u 


d 


s 


c 


h 


u 


222.8 


P° 

768(768) 
714(135) 


P+ 
773(768) 


K*+ 
910(892) 


ir u 

2110(2007) 


B* + 
5712(5325) 


d 


236.2 


658(140) 


to 

782(782) 

668(549) 
V 


914(896) 


D*~ 
2109(2010) 


5* u 
5709(5325) 


s 


427.2 


825(494) 
K~ 


831(498) 
K° 


1019(1019) 

953(958) 
V 


D*- 
2156(2110) 


B? 
5735( — ) 


c 


1701.3 


2079(1865) 
D° 


2078(1869) 
D+ 


2131(1969) 
Df 


3097(3097) 
3082(2980) 

Vc 


B* + 
6502( — ) 


b 


5328.2 


5701(5278) 
B~ 


5698(5279) 
B° 


5726(5375) 


6495( — ) 
Be 


T 

9460(9460) 

9455( — ) 
Vb 
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with the experimental values to the extent the latter are known. The present model predicts for 
example 

M(B±) = 6495 MeV , M(B* ± ) = 6502 MeV . (58) 

By and large the agreement is remarkably good. The heavy meson masses are reproduced quite 
well, but the agreement is not quantitative everywhere, particularly not for those hadrons with one 
light quark (u or d). In judging this agreement one should keep in mind, (1) that such a table, 
in which all hadrons have been calculated from one and the same model, has hitherto not been 
prepared; and (2) that the light mesons like the quasi-pions should not be calculated by a crude 
potential model like the retarded Schrodinger equation. The smallness condition actually tells us 
that these hadrons probably are systems in which the constituents move highly relativistically. 
Thus far, we have at hand no simple paradigms for such a kinematic situation. Solving directly the 
momentum-space integral equations might therefore be the only way. 

All in all, with all due respect to the work with potential models and with lattice gauge theory, 
the agreement between the empirical data and the present first attempt to relate them on trial and 
error to an effective, QCD-inspired Hamiltonian is in fact not so much worse; particularly in view 
of the absence of any free parameters and the simplicity of the approach. No doubt, the various 
simplifications can be improved in the future. 



6 Summary and Conclusions 

The full many-body front-form Hamiltonian, evaluated for QCD in the light-cone gauge A + = 0, 



had been reduced in preceeding work [21] to a manifestly gauge invariant effective Hamiltonian 
which acts only in the space of one quark and one antiquark. Particularly, no Tamm-Dancoff 
Fock-space truncations had to be made, nor was it necessary to rely on perturbation theory by 
assuming a small coupling constant. The present work is motivated by the question why and how 
such a simple structure like the resulting integral equation in the transversal momenta and the 
longitudinal momentum fraction can account for the complexities of hadronic phenomenology. In 
particular we have wondered to which extent one can understand the masses of the pseudoscalar 
and vector mesons with no other input than the flavor quark masses of the constituent quarks. 



In this first study of such a structure, which actually was preceeding |27|, |2C|] the more rigorous 
derivation [21], we replace the running coupling constant, which absorbes the many-body ampli- 



tudes of the full theory in a well defined way, by the suitably adjusted phenomenological version 
of Richardson ]22]. At the least, the latter interpolates smoothly between asymptotic freedom and 
infrared slavery. Its only free parameter is fixed by a fit to the strong coupling constant at the Z 
mass. 

For the future, we have in mind mainly two improvements: (1) the explicit calculation of the 
running coupling constant by a renormalization group analysis, and (2) an explicit solution of the 
integral equation in light-cone variables Eq.([|). It should be applied to mesons whose constituent 
quarks havevery different masses, such that the structure functions including the contributions from 
higher Fock states can be calculated from a covariant theory. This could be done in such a way 
that the relation to the existing phenomenological models can be seen explicitly. 

With these future applications in mind, we have not hesitated to perform in Section 3 a number 
of basically trivial and straightforward calculation and to transcribe the front-form integral equation 
into the intuitively easier accessible form of usual momenta. The major impact of very different 
quark masses is then absorbed into the familiar reduced mass, and all integration variables have 
the same domain of validity. This is not unimportant for the practitioner who actually wants 
to get out numbers from his theory. This virtue does not seem to be common place anymore, 
unfortunately. As a wonderful and not intended side effect, it turns out that the rotationally 
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only covariant equation on the light-cone can be transformed into a rotationally invariant integral 
equation Eq.(|36|) in usual momentum space. All factors which seem to violate strict rotational 
invariance can be absorbed into the wavefunction, Eq.(^). One looks forward to see numerical 
solutions to these equations. 

But in our aim to relate the basically exact formalism with its connection to Lagrangian QCD 
to the usual configuration space where our intuition is at home, we went a step further and tried to 
Fourier transform the integral equations. We have been unable to do this, by formal mathematical 
reasons. Rather we had to discourse to an approximation to which we refer to as semi-relativistic. 
The resulting retarded Schrodinger equation Eq.(^Tj) has the amazing property of looking like 
a conventional Schrodinger equation with velocity-dependent interactions and still being a fully 
covariant equation. It should be obvious that the transition from the front form to the usual instant 
form momenta and the subsequent Fourier transform to configuration space does not change the 
basic feature of the light-cone Hamiltonian to be manifestly frame independent. Would one be able 
to perform the necessary Fourier transforms in closed form, this statement could be phrased even 
more rigorously. 

One should emphasize that the retarded Schrodinger equation has no free parameter, since 
coupling constant and quark masses have to be determined from the experiment. Fitting the 5 
quark flavor masses to the 5 qq-vector mesons exhausts all freedom. The rest is structure: The 
25 remaining pseudoscalar and vector masses are then predicted and presented in Table ||. In 
comparison with the experimental data, they are not much worse than those from conventional 
phenomenological models |]2^ , 30 1, or from heavy quark symmetry [34], or even from lattice gauge 
calculations |35|, 37 1, in particular when keeping in mind the very rough and simple approxi- 
mations applied. The pions are reproduced more than poorly and remain mysterious particles like 
in every other model not specially designed for them. The mesons with one light quark do not 
yet meet the tough standard of the phenomenological models. The latter two aspects are possibly 
related to each other. 

Conclusion: If such a poor model can do so well one must be on the right track. It seems that 
the front form Hamiltonian approach applied to quantum chromodynamics has made a big step 
forward. Intensified efforts are justified. 
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